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Abstract: It is a common impression that by only setting the maximum occupation 
number to infinity, which is the demand of the indistinguishability of bosons, one can 
achieve the statistical distribution that bosons obey — the Bose-Einstein distribution. 
In this letter, however, we show that only with an infinite maximum occupation number 
one cannot uniquely achieve the Bose-Einstein distribution, since in the derivation of the 
Bose-Einstein distribution, the problem of iterated limit is encountered. For achieving 
the Bose-Einstein distribution, one needs to take both the maximum occupation number 
and the total number of particles to infinities, and, then, the problem of the order of 
taking limits arises. Different orders of the limit operations will lead to different statistical 
distributions. For achieving the Bose-Einstein distribution, besides setting the maximum 
occupation number, we also need to state the order of the limit operations. 
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1 Introduction 

The indistinguishability of identical particles demands that for Bose-Einstein statistics, a 
quantum state can be occupied by any number of particles, and for Fermi-Dirac statistics, a 
quantum state can be occupied by only one particle. Starting from this point, as a common 
impression, in statistical mechanics, one can obtain statistical distributions by only setting 
the various values of the maximum occupation number n: by n — > oo, one achieves the 
Bose-Einstein distribution, and by n = 1, one achieves the Fermi-Dirac distribution. In 
this letter, however, we show that, contrary to this common impression, only from a 
maximum occupation number, one cannot uniquely construct the statistical distribution 
in the Bose-Einstein case, though such a treatment works well in the Fermi-Dirac case. 
This difference between Bose-Einstein and Fermi-Dirac cases comes from the fact that the 
maximum occupation number in the Bose-Einstein case is oo, but in the Fermi-Dirac case 
is a finite number. 

In section [21 we point out that there exists not only one possibility to take limits in the 
derivation of the Bose-Einstein distribution. In sections [3] and [U we discuss the various 
orders of taking limits and show which one will lead to the Bose-Einstein distribution. 
The conclusion is summarized in Sec. [5j 
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2 The order of taking limits 



The reason why the Bose-Einstein distribution cannot be uniquely achieved by only set- 
ting n — > oo is that there are two infinities involved in the derivation of the Bose-Einstein 
distribution — the maximum occupation number and the total number of particles. The 
infinite maximum occupation number is the demand of the indistinguishability of bosons; 
the infinite total number of particles, i.e., the thermodynamic limit, guarantees the exten- 
sivity of a thermodynamic system. As a result, during the derivation, we need to perform 
two limits, lim and lim , where (N) denotes the mean total number of particles in the 

grand canonical ensemble. This means that in the Bose-Einstein case, we have to first 
make a choice of the order of the limit operations. 

For simultaneously taking the maximum occupation number n and the total number 
of particles (N) into account, we start with Gentile statistics [H[2] in which the maximum 
occupation number n can take on any value. Different values of n correspond to different 
kinds of statistics, and the Bose-Einstein case is a special case of Gentile statistics with 
n — > oo [3J. The Gentile distribution reads 

f 1 " + 1 m 

z -l e e/(*T)_l [ z -l eE /(*T)] n + 1 _l' KJ 

where T is the temperature, e is the energy, and z is the fugacity. In the grand canonical 
ensemble, the fugacity z is determined by [3] 



where A = h/^/2-KmkT is the mean thermal wavelength, 

h/2 (z) = 53/2 (*) - J— fl 3 / 2 {z n+l ) , 
V n + 1 

and $3/2 (z) is the Bose-Einstein integral. In principle, from eq. ([2]), one can solve the 
expression of z as 

z = w(T,n,(N)), (3) 
and then write the distribution function as 

f G (T,n,(N)) 



w(T,n, (N))" 1 e £ /(* T ) 
n + 1 



1 n+l 

w(T,n,(N)y 1 e £ /( kT ) -1 



(4) 



There are two infinities, n and (N), in the distribution function fa (T,n, (N)). For 
recovering the Bose-Einstein distribution, we need to take both n and (N) tend to oo in 
eq. dH). However, there are three possible ways to take the limits: two iterated limits 

lim lim f G (T,n,{N}) (5) 

(TV)— >oo n-^oo 
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and 



lim lim f G (T,n,(N)), (6) 

n-*oo (N)-*oo 



and the double limit 

lim fa (T, n, (N)) . (7) 

n—>oo 
{N)->oc 

The first thing is thus to choose the way how to take the limit from the above three. Nev- 
ertheless, the indistinguishability of bosons only requires an infinite maximum occupation 
number, but does not say anything about the order of the limits. Generally speaking, dif- 
ferent orders of the two limit operations will lead to different results [HE]. At first sight, 
one may think that a natural choice is the double limit (|7|). However, in the following, we 
will show that the double limit ([7]) does not exist, since the iterated limits (0) and ([U]) exist 
but are not equal to each other. It is the iterated limit ([5]) that leads to the Bose-Einstein 
distribution. That is to say, for obtaining the Bose-Einstein distribution, besides setting 
n — > oo, the demand of indistinguishability, we also need an additional condition that the 
order of the iterated limit should be lim lim . 

3 lim lim fc' the Bose-Einstein distribution 

We first discuss the case of the iterated limit lim lim fc- In this case, the limit lim 
is first performed. 

In principle, for obtaining the statistical distribution under the limit of n — > oo, we 
need to calculate the explicit expression for the fugacity z, substitute it into eq. ([T]) to 
achieve eq. (HJ), and then perform the limit lim in eq. (jH). However, it is difficult to 

n— >oo 

obtain the expression of z, so we will do this in an indirect way. 

We first determine the maximum value of the fugacity z at n — > oo. Since — 
z -(n+i) — f is a monotonically increasing function of z, from eq. ([2]), the maximum value of 
z appears at T = 0, or, A — > oo, i.e., 



n^oo 



(AT) = lim 

The solution of eq. (JH]) is 



ra + l 



r -l _ i -(n+1) 

Zmax J- z max 



(8) 



C) 



— (AT) + 1 

Based on this result, we can take the limit lim on eq. (J2|): 

n— >oo 

(N) = ^g 3/2 (z) + —}— [ , (10) 

or 

1 = -^93/2 (*) + Tj^rp (ii) 
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where p = (N) /V is the particle number density. 

Now, we can take the limit lim and obtain the range of value of the fugacity 

{N)-+oo 

0<z<l. (12) 

This is just the Bose-Einstein case [6]. The corresponding statistical distribution is the 
Bose-Einstein distribution: ^ 

f = z -l e e/(kT) _i ( 13 ) 

The above result shows that if the order of the limits is chosen as lim lim fa, one 

(A 7 )— >oo n-^oo 

achieves the Bose-Einstein distribution. 



4 lim lim fc- beyond the Bose-Einstein distribution 

n->oo (A/)— kx) 

An alternative choice for the iterated limit is lim lim , i.e., the limit lim is first 
performed. 



Rewrite eq. ([2j) as 

1 1 h ( \ -L 1 ' 1 " + ' 



1 z -(n+l) _ ! 



(14) 



Here — is the number of the particles occupying the ground state, whose 

maximum value is the maximum occupation number n. When taking the limit lim , 

{N)-*oo 

the ground-state term vanishes. Consequently, taking the limit lim on eq. (1141) gives 

p\ 3 = hy 2 (z) . (15) 

(z) is a monotonically increasing function of z and < ft.3/2 (z) < 00 [3j, so the range 
of the value of z is < z < 00. 

Next, take the limit lim on eq. (|15|) , 

n^oo 

After taking the limit lim , eq. (|14p becomes eq. ([15p . The left-hand side of eq. 

(|15[) can take on any value from to 00, so the fugacity z can take on the value that is 
greater than 1 due to the fact that h s / 2 (z) is a monotonically increasing function of z and 
when z = 1, lim /13/2 (z) = £ (3/2), where £ (3/2) is the zeta function. 

n— >oo ' 

The distribution corresponding to 2 > 1 must not be the Bose-Einstein distribution 
because in the Bose-Einstein case, < z < 1. Concretely, taking lim in the distribution 

n— +00 

function gives 

f G =\ *-VA fcT )-l' e> / i ( T )' / 16 s 
\ 00. e < p{T) , 

where ^ is the chemical potential. Notice that before taking lim , the limit lim has been 

ra^oo (A") — >oo 

already taken. Eq. (|16p is a statistical distribution with an infinite maximum occupation 
number, i.e., n — > 00, but is not the Bose-Einstein distribution. 
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The above results show that both the iterated limits lim lim f@ and lim lim fc 
exist, but are not equal to each other. As a result, the double limit lim Jq does not 

n—>oo 
(N)^oo 

exist. 

5 Conclusions 

In conclusion, in the derivation of the Bose-Einstein distribution, two infinities are en- 
countered. Thus, we have to choose the way how to take the limits. Different orders 
of the limit operations give different statistical distributions. In the common derivation 
of Bose-Einstein distribution, the order of the iterated limit has been unstatedly chosen 
as lim lim . That is to say, to achieve the Bose-Einstein distribution, besides setting 

(AT)— s-oo n->oo 

the maximum occupation number to infinity, one also needs to choose the order of taking 
limits. 
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